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Abstract 

This is the first paper of a series devoted to the study of the cosmological imphca- 
tions of the parallel mirror world with the same microphysics as the ordinary one, but 
having smaller temperature, with a limit set by the BBN constraints. The difference 
in temperature of the ordinary and mirror sectors generates shifts in the key epochs for 
structure formation, which proceeds in the mirror sector under different conditions. 
We consider adiabatic scalar primordial perturbations as an input and analyze the 
trends of all the relevant scales for structure formation (Jeans length and mass. Silk 
scale, horizon scale) for both ordinary and mirror sectors, comparing them with the 
CDM case. These scales are functions of the fundamental parameters of the theory 
(the temperature of the mirror plasma and the amount of mirror baryonic matter), 
and in particular they are influenced by the difference between the cosmological key 
epochs in the two sectors. Then we used a numerical code to compute the evolution 
in linear regime of density perturbations for all the components of a Mirror Universe: 
ordinary baryons and photons, mirror baryons and photons, and possibly cold dark 
matter. We analyzed the evolution of the perturbations for different values of mirror 
temperature and baryonic density, and obtained that for x — T' /T less than a typical 
value Xeq, for which the mirror baryon-photon decoupling happens before the matter- 
radiation equality, mirror baryons are equivalent to the CDM for the linear structure 
formation process. Indeed, the smaller the value of x, the closer mirror dark matter 
resembles standard cold dark matter during the linear regime. 



1 Introduction 



The present cosmological observations strongly support the main predictions of the infla- 
tionary scenario: first, the Universe is flat, with the total energy density very close to the 
critical fig ~ 1, and second, primordial density perturbations have nearly flat spectrum, 
with the scalar spectral index n ~ 1. The non-relativistic matter gives only a small frac- 
tion of the present energy density, about Q.m ^ 0.25, while the rest is attributed to the 
vacuum energy (cosmological term or dark energy) VL\ ~ 0.75. The fact that Vtm and Vt\ 
are of the same order gives rise to the so called cosmological coincidence problem: why we 
live in an epoch when pm ~ Ph-, if in the early Universe one had pm ^ PK and in the late 
Universe one would expect pm ^ PA? The answer can be only related to an anthropic 
principle: the matter and vacuum energy densities scale differently with the expansion of 
the Universe, pm oc and p\ oc const., hence they have to coincide at some moment, 
and we are just happy to be here. Moreover, for substantially larger no galaxies could 
be formed and thus there would not be anyone to ask this question. 

On the other hand, the matter in the Universe has two components, visible and dark: 
ilm = ^fe + ^d- The visible matter consists of baryons with Of, ~ 0.04, while the dark 
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matter with ~ 0.2 is constituted by some hypothetic particle species very weakly 
interacting with the observable matter. It is a tantalizing question why the visible and 
dark components have so close energy densities? Clearly, the ratio Pd/ Pb does not depend 
on time as far as with the expansion of the Universe both and pd scale as oc a~^. 

In view of the standard cosmological paradigm, there is no good reason for having 
~ ilfi, as far as the visible and dark components have different origins. The density of 
the visible matter is pb = M^nt, where — 1 GeV is the nucleon mass, and is the 
baryon number density of the Universe. The latter should be produced in a very early 
Universe by some baryogenesis mechanism, which is presumably related to some B and 
CP-violating physics at very high energies. The baryon number per photon rj = rib/n^ is 
very small. Observational data on the primordial abundances of light elements and the 
recent results on the cosmic microwave background radiation (CMB) anisotropics nicely 
converge to the value ~ 6 x 10"^''. 

As for dark matter, it is presumably constituted by some cold relics with mass M and 
number density n^, and pd = Mud- The most popular candidate for cold dark matter 
(CDM) is provided by the lightest supersymmetric particle (LSP) with Misp ~ 1 TeV, 
and its number density nisp is fixed by its annihilation cross-section. Hence pf, ~ Plsp 
requires that nb/riLSP ~ Mlsp/Mn and the origin of such a conspiracy between four 
independent parameters is unclear: the value M^r is fixed by the QCD scale while Mlsp 
is related to the supersymmetry breaking scale, rn, is determined by B and CP violating 
properties of the particle theory at very high energies whereas nisp strongly depends on 
the supersymmetry breaking details. Within the parameter space of the MSSM (Minimal 
Supersymmetric Standard Model) it could vary within several orders of magnitude, and 
moreover, in either case it has nothing to do with the B and CP violating effects. The 
situation looks even more obscure if the dark component is related e.g. to the primordial 
oscillations of a classic axion field, in which case the dark matter particles constituted by 
axions are superlight, with mass <C 1 eV, but they have a condensate with enormously 
high number density. 

A new sight on the problem of dark matter can be given by the concept of mirror 
world. This idea was studied in ref. and then its phenomenological and cosmological 
implications discussed in several later papers. In particular, after its first applications 
to non-baryonic dark matter the mirror matter hypothesis has attracted a significant 
interest over the last years and has been invoked in many physical and astrophysical ques- 
tions: large scale structure of the Universe galactic halo 5 , MACHOs gamma 
ray bursts [7] , flavour and CP violation |H] orthopositronium lifetime (HI j interpreta- 
tion of dark matter detection experiments , meteoritic event anomalies |12| , close-in 
extrasolar planets I14j . Pioneer spacecraft anomalies jl5| . 

Mirror world can be considered as a hidden sector of particles and interactions which 
are exactly the same as in our visible world^, and coupling to our world essentially by 
gravity. Besides the gravity, the two sectors could communicate by other means. In par- 
ticular, ordinary photons could have kinetic mixing with mirror photons ^7j, or ordinary 
(active) neutrinos could mix with mirror (sterile) neutrinos |18j . 

If the mirror sector exists, then the Universe along with the ordinary (O) photons, 
neutrinos, baryons, etc. should contain their mirror (M) partners. One could naively 
think that due to mirror parity the ordinary and mirror particles should have the same 
cosmological abundances and hence the O and M sectors should have the same cosmological 

^The mirror parity could be spontaneously broken and the weak interaction scales ((/>) = v and {(f)') — v' 
could be different, which leads to somewhat different particle physics in the mirror sector In this paper 
we treat only the simplest case, v — v' , in which the mirror sector has exactly the same physics as the 
ordinary one. 
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evolution. However, this would be in conflict with the Big Bang nucleosynthesis (BBN) 
bounds on the effective number of extra light neutrinos, since the mirror photons, electrons 
and neutrinos would give a contribution to the Hubble expansion rate equivalent to AA'^^ ~ 
6.14 Therefore, even if their microphysics is the same, initial conditions are not 

the same in both sectors, and in the early Universe the M system should have a lower 
temperature than ordinary particles. This situation is plausible if the two systems are 
born with different temperatures (the post-inflationary reheating temperature in the M 
sector is lower than in the visible one jl9l l2Uj ) . they interact very weakly, so that they 
do not come into thermal equilibrium with each other (condition automatically fulfilled if 
the two worlds communicate only via gravity) , and both of them expand adiabatically. So 
they evolve independently and their temperatures remain different at later stages. 

At present, the temperature of ordinary relic photons is T ~ 2.72 K, and the mass 
density of ordinary baryons constitutes about 5% of the critical density. Mirror photons 
should have smaller temperature T' < T, so their number density is n'^ = x^n^, where 
X = T' /T. ^ This ratio is a key parameter in our further considerations since it remains 
nearly invariant during the expansion of the Universe. The BBN limit on AN^, implies 
the upper bound x < O.GA ANu^'^ . If we assume AN^, ^ 1, we obtain x ^ 0.64 as a limit 
for the temperature of the mirror sector |21j|. The dependence of x on ANi, is indeed very 
weak, so even a value x ^ 0.7 could be compatible with the current experimental status. 
As for mirror baryons, their number density can be larger than n^,, and if the ratio 
P = nyub is about 5 or so, they could constitute the dark matter of the Universe. 

In this view, the mirror world idea could give a new twist to dark matter problem. 
Once the visible matter is built up by ordinary baryons, then the mirror baryons could 
constitute dark matter in a natural way. They interact with mirror photons, however they 
are dark in terms of the ordinary photons. The mass of M baryons is the same as the 
ordinary one, = Mat, and so we have /? = QyClb- In addition, as far as the two sectors 
have the same particle physics, it is natural to think that the M baryon number density 
n[ is determined by the baryogenesis mechanism which is similar to the one which fixes 
the O baryon density rib. Thus, one could question whether the ratio (3 could be naturally 
order 1 or somewhat bigger. 

There are several baryogenesis mechanisms as are GUT baryogenesis, leptogenesis, 
electroweak baryogenesis, etc. At present it is not possible to say definitely which of 
the known mechanisms is responsible for the observed baryon asymmetry in the ordinary 
world. However, it is most likely that the baryon asymmetry in the mirror world is 
produced by the same mechanism and moreover, the properties of the B and CP violation 
processes are parametrically the same in both cases. But the mirror sector has a lower 
temperature than ordinary one, and so at epochs relevant for baryogenesis the out-of- 
equilibrium conditions should be easier fulfilled for the M sector. 

In particular, we know that in certain baryogenesis scenarios the M world gets a larger 
baryon asymmetry than the O sector, and it is pretty plausible that f3 > 1 ^21 . This 
situation emerges in a particularly appealing way in the leptogenesis scenario due to the 
lepton number leaking from the O to the M sector which leads to nj, > n^, and can thus 
explain the near coincidence of visible and dark components in a rather natural way |22j . 

Given its consistency with all the experimental situation so far, it's important to 
construct a complete theory of cosmology with mirror dark matter and test this scenario. 
Thus the aim of this work (described in this paper and the next one of this series, hereafter 
referred to as Paper 2 |23]) is contained in a question: "is mirror matter still a reliable dark 
matter candidate?" Its emergence arises also from the astrophysical problems encountered 

■^From now on all quantities of the mirror sector will be marked by ' to distinguish from the ones 
belonging to the observable or ordinary world. 
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by the standard candidate, the cold dark matter (CDM), as for example the central galactic 
density profiles or the number of small satellites. 

In this first paper, we focus on the structure formation theory in linear regime, ana- 
lyzing the trends of scales in both sectors and comparing them with the CDM case, and 
finally showing the evolution of perturbations in all the components of a Mirror Universe^ 
(ordinary and mirror photons and baryons, and possibly CDM). In doing so, we largely 
extend our previous studies of ref. 0], analyzing the details of the structure formation 
and evolution in a mirror scenario. The paper is organized as follows. In section 2 we 
show that, due to the temperature difference, in mirror sector the key epochs for structure 
formation occur at different redshifts than in the observable sector. The section 3 ana- 
lyzes the relevant scales for structure formation (sound speed, Jeans length and mass. Silk 
mass) as functions of the temperature and the baryonic density of the mirror sector. The 
main differences with respect to ordinary baryonic and CDM scenarios are also discussed. 
In section 4 we study the possible mirror scenarios for the growth of primordial pertur- 
bations. Section 5 is devoted to the study of the temporal evolution of perturbations in 
all the ordinary and mirror components. Finally, in section 6 we briefly summarize our 
findings. 

2 Mirror baryons as dark matter 

Since it interacts only gravitationally with our ordinary sector, the mirror matter is a nat- 
ural dark matter candidate. At present there are observational evidences that dark matter 
exists and it's density is about 5 times that of ordinary baryons, but in previous section 
we said that this is not a problem for mirror scenario, because the same microphysics does 
not imply the same initial conditions in both O and M sectors. 

In the most general context, the present energy density contains relativistic (radiation) 
component Qr, non-relativistic (matter) component 0^ and the vacuum energy density 
Oa. According to the inflationary paradigm the Universe should be almost flat, f^o = 
Qrri + + ~ 1; which well agrees with the results on the CMB anisotropy. 

If we consider now a Mirror Universe, i.e. a Universe made of two sectors, ilr and 
represent the total amount of radiation and matter of both ordinary and mirror sectors: 
Q,r = {^r)o + {^r)M and Qm = {^m)o + {^m)M- The two parameters describing the 
mirror sector are the ratio of the temperatures of two sectors x and the relative amount 
of mirror baryons compared to the ordinary ones /?: 

x = ^< 0.64 ^ = § ^ ^ • 

In the above expressions, the first limit comes from the BBN bound on the effective number 
of extra light neutrinos (see previous section) and the second one from the hypothesis that 
mirror baryonic contribution to dark matter has cosmological relevance. In the context 
of our model, as explained in |2J, the relativistic fraction is represented by the ordinary 
and mirror photons and neutrinos, and, using the expression for the ordinary degrees of 
freedom in a Mirror Universe, g(T) = g^(T){l + x'^), and the value of the observable 
radiation energy density {^r)o — 4.2x10^^, it is given by 

VLr = 4.2 X 10"^ /i"^ (1 + x^) ~ 4.2 x 10"^ /i"^ , (2) 

^The expression "Mirror Universe" is clearly misleading, since one could think that there is another 
Universe, while it is just one, but made of two sectors. Nevertheless, this expression is sometimes used to 
shortly refer to this scenario. 
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where the contribution of the mirror species, expressed by the additional term x"^, is 
neghgible in view of the BBN constraint x < 0.64. As for the non-relativistic component, 
it contains the O baryon fraction Qf, and the M baryon fraction 0,'^ = f3Qh, while the 
other types of dark matter could also present. Obviously, since mirror parity doubles 
all the ordinary particles, even if they are "dark" (i.e., we are not able to detect them 
now), whatever the form of dark matter made by some exotic ordinary particles, there 
will exist a mirror partner made by the mirror counterpart of these particles. In the 
context of supersymmetry, the CDM component could exist in the form of the lightest 
supersymmetric particle (LSP). It is interesting to remark that the mass fractions of the 
ordinary and mirror LSP are related as ^'isp — ^^LSP- In addition, a significant HDM 
component could be due to neutrinos with order eV mass. The contribution of the 
mirror massive neutrinos scales as ^l'^ = x^Oi, and thus it is irrelevant. In any case, 
considering the only CDM component, which is now the preferred candidate, we can 
combine both the ordinary and mirror components, since their physical effects are exactly 
the same. Thus we consider a matter composition of the Universe expressed in general by 

nrn = ^b + ^'b + ^CDM = ^b{l + P) + ^CDM ■ (3) 

The important moments for the structure formation are related to the matter-radiation 
equality (MRE) epoch and to the plasma recombination and matter-radiation decoupling 
(MRD) epochs in both sectors. The MRE occurs at the redshift 

which is always smaller than the value obtained for an ordinary Universe, but approximates 
it for low x (see fig. If we consider only ordinary and mirror baryons and photons, we 
find 

Pb{l + P) , . ^ P'b{l + P- '' 

p^{i + x^y p'^{i + x~ , 

where the baryon and photon densities refer to the redshift z. This implies that, with the 
addition of a mirror sector, the matter-radiation equality epoch shifts toward earlier times 



l + -e.= ^^^)^(l + .) = ^n^^(l + .), (5) 
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as 

(l + x^ 

The MRD, instead, takes place in every sector only after the most of electrons and 
protons recombine into neutral hydrogen and the free electron number density dimin- 
ishes, so that the interaction rate of the photons = n^aT = Xgrju^aT drops below 
the Hubble expansion rate H{T), where cjt is the Thomson cross section, X^, = ne/ni, is 
the fractional ionization, and rj = n^/n^ is the baryon to photon ratio. In condition of 
chemical equilibrium, X^. is given by the Saha equation, which for X^. <^ 1 reads 

^e-(l-n)^/^^(-)"'^%-^/^^ (7) 

r/^/^ \me/ 

where B = 13.6 eV is the hydrogen binding energy and I4 is the ^He abundance. Thus 
we obtain the familiar result that in our Universe the MRD takes place in the matter 
domination period, at the temperature T^ec — 0.26 eV, which corresponds to redshift 

1 + 2dec = 7dec/To ~ 1100. 

*In the most general case, where there is also some other dark matter component, as the CDM, the 
parameter /3 could be the sum of two terms, /3 + /3cdm- 
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The MRD temperature in the M sector T^^^ can be calculated following the same lines 
as in the ordinary one. Due to the fact that in either case the photon decoupling occurs 
when the exponential factor in eq. Q becomes very small, we have T^^^ ~ Tjeo up to small 
corrections related to ij' , I4 different from tj, ^4 (for more details see ref. Hence, 
considering that T' = x ■ T , we obtain 

1 + ~ X-\l + Zdec) ^ 1.1 • lO^X-l , (8) 

SO that the MRD in the M sector occurs earher than in the ordinary one. Moreover, 
comparing eqs. (jlj) and 0, which have different dependences on x, we find that, for x 
smaller than a typical value Xeq expressed by 

Xeq « 0.046(J^„/l2)-l , (9) 

the mirror photons would decouple yet during the radiation dominated period (see fig.^. 
Assuming, e.g., the value l^m^^ = 0.135 (the result of a recent WMAP fit we obtain 
Xeq ~ 0.34, which indicates that below about this value the mirror decoupling happens in 
the radiation dominated period, with consequences on structure formation (as we will see 
in the following sections and in paper 2 j23j ) . 




0.1 1.0 



Figure 1: The M photon decoupling redshift 1 + as a function of x (thick solid). The horizontal 
thin solid line marks the ordinary photon decoupling redshift 1 + Zdcc = 1100. We also show the 
matter-radiation equality redshift l+Zcq (dash) and the mirror Jeans- horizon mass equality redshift 
1 + z'^ (dot) for the case ilmh^ — 0.135. The shaded area x > 0.64 is excluded by the BBN limits. 

We have shown that mirror baryons could provide a significant contribution to the 
energy density of the Universe and thus they could constitute a relevant component of 
dark matter. Immediate question arises: how the mirror baryonic dark matter (MBDM) 
behaves and what are the differences from the more familiar dark matter candidates as 
the cold dark matter (CDM), the hot dark matter (HDM), etc. In next sections we discuss 
the problem of the cosmological structure formation in the presence of M baryons as a 
dark matter component. 

3 Mirror baryonic structure formation 

In this section, we extend the linear structure formation theory (for the standard scenario, 
see refs. and j26j ) to the case of dark matter with a non-negligible mirror baryonic 
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component. 

In a Mirror Universe we assume that a mirror sector is present, so that the matter 
is made of ordinary baryons (the only certain component), non-baryonic (dark) matter, 
and mirror baryons. Thus, it is necessary to study the structure formation in all these 
three components. We proceed in this way: first of all, we recall the situation for ordinary 
baryons, and then we compute the same quantities for mirror baryons, comparing them 
with each other and with the CDM case. 

In general, when dealing with the pre-recombination plasma, we distinguish between 
two types of perturbations, namely between "isoentropic" (adiabatic) and "entropic" 
{isocurvature or isothermal) modes |27| . while after matter-radiation decoupling pertur- 
bations evolve in the same way regardless of their original nature. 

In this paper we study only adiabatic perturbations, which are today the preferred 
perturbation modes, and we leave out the isocurvature modes, which could also have a 
contribution, but certainly cannot be the dominant component [2H1- Here we remember 
only that an adiabatic perturbation satisfies the condition for adiabaticity 

drr. = ^Sr , (10) 

which, defining, as usual, 5 = {5p/p), relates perturbations in matter {6m) and radiation 
{6r) components. 

We will now consider cosmological models where baryons, ordinary or mirror, are the 
dominant form of matter. Thus, it is crucial to study the interaction between baryonic 
matter and radiation during the plasma epoch in both sectors, and the simplest way of 
doing it is by looking at models containing only these two matter components. 

According to the Jeans theory [29^ , the relevant scale for the gravitational instabilities 
is characterized by the Jeans scale (length and mass), which now needs to be defined in 
both the ordinary and mirror sectors. Then, we define the ordinary and mirror Jeans 
lengths as 

V L^Pdom V t^Pdom 

where pdom is the density of the dominant species, and Vg, v!, are the sound speeds, and 
the Jeans masses as 

Mj = l^p, = ^p, (X.,)' = Inp', (I)' , (12) 

where the density is now that of the perturbed component (ordinary or mirror baryons). 

3.1 Evolution of the adiabatic sound speed 

Looking at the expressions of the Jeans length (|llj) . it is clear that the key issues are 
the evolutions of the sound speeds in both sectors, since they determine the scales of 
gravitational instabilities. Using the definition of adiabatic sound speed, we obtain for the 
two sectors 

where w and w' are relative respectively to the ordinary and mirror equations of state 
p = wp and p' = w' p' . 

First of all, we consider the standard case of a Universe made of only one sector. In a 
mixture of radiation and baryonic matter the total density and pressure are p = p^ + ph 
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and p c:i p-y 
given by 



p-y/3 respectively (recall that pb 



dp 
dp 



1/2 



0). Hence, the adiabatic sound speed is 

1/2 



(14) 



where we have used the adiabatic condition H1U|) . In particular, using the scaling laws 
Pm oc /Oom(l + z)^ and py oc ^07(1 + z)^, together with the definition of matter-radiation 
equality (where we now consider only baryons and photons), we obtain 



Vsiz) ~ 



1 

7! 



1 + 



3 /l + z, 



eq 



l + Z 



1-1/2 



(15) 



In fact, the relation above is valid only for an ordinary Universe, and it is an approximation, 
for small values of x and the mirror baryon density (remember that (3 = $7^/17;,), of the 
more general equation for a Universe made of two sectors of baryons and photons, obtained 
using eqs. ((T^ and © and given by 



Vs{z) ~ 



1 
73 



1 + 



3 / 1 + 



'l + z, 



eq 



l + Z 



-1/2 



(16) 



In the most general case, the matter is made not only of ordinary and mirror baryons, but 
also of some other form of dark matter, so the factor 1 + /3 is replaced by 1 + (3 + Pdm , 
where Pdm = {^m — — ^'b) l^b- The presence of the term [(H-x^)/(l + /3)] in equation 
(|16|) is linked to the shift of matter-radiation equality epoch @, thus it only balance this 
effect, without changing the value of the sound speed computed using eq. (|15|) . 

The mirror plasma contains more baryons and less photons than the ordinary one, 
I3pb and p'^ = x'^p^. Then, using eqs. (|T^ and ©, we have 



P'b 



iz) 




-1/2 



1 
73 




1-1/2 



(17) 



Let us consider for simplicity the case when dark matter of the Universe is entirely due 
to M baryons, — ^'b (i-6., P ^ 1). Hence, for the redshifts of cosmological relevance, 
z ~ Zeqi we have vi. ~ 2x^/3, which is always less than Vg ~ l/\/3 (some example: if 
X = 0.7, vi. « 0.5 • Vg] if X = 0.3, vi. ^ 0.1 • v^). In expression ((T7)) it is crucial the presence 
of the factor [(1 + x~^)/(l + which is always greater than 1 (given the bounds of 

eqs. on the parameters), so that vi. < Vg during all the history of the Universe, and 
only in the limit of very low scale factors, a <C Oeq, we obtain t>g ~ ~ l/\/3- As we will 
see in the following, this has important consequences on structure formation scales. 
Now we define Qiy^y Si's the scale factor corresponding to the redshift 



(1 + Zb-y) 



{nb/n^) = 3.9 • lo^ifibh^ 



(18) 



Since 1 + z^ec — 1100, ordinary baryon-photon equipartition occurs before recombination 
onlyif Qft/i^ > 0.026 (which seems unlikely, given its current estimates). According to what 
found in § 121 in the mirror sector the scale of baryon-photon equality a[,^ is dependent on 
X and it transforms as 



/ 

Ob-) 



(19) 



If we remember the definition of the quantity Xeq ~ 0.046(17^^ ) , we have that for x > 
the decoupling occurs after equipartition (as in the ordinary sector for ri^/i^ > 0.026), 



X, 



while for x < Xcq it occurs before (as for l^^/i^ < 0.026). 
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Figure 2: The trends of the mirror sound speed (vi.) as a function of the scale factor, compared 
with the ordinary sound speed {vs). The ordinary model has flhh^ — 0.08, while the mirror model 
has X = 0.6 and (3 = 2. Are also reported all the key epochs: photon-baryon equipartition and 
decoupling in both sectors, and the matter-radiation equality. 

Regardless of which sector we are considering, in the radiation era ^ ph, ensuring 
that Vs — l/\/3. In the interval between equipartition and decoupling, when ^> p^, 
eq. ((TH) gives ~ Ap^/Spi, oc a^'^^'^. After decoupling there is no more pressure equi- 
librium between baryons and photons, and Vg is just the velocity dispersion of a gas of 
hydrogen and helium, Vs oc a~^. If Ofc/i^ < 0.026 or x < Xeq (according to what sector 
we consider), photon-baryon equipartition occurs after decoupling, and the intermediate 
situation does not arise. 

It follows that, by taking care to interchange with a'^^ and adec with a^ec have 
for the sound speed the same trends with the scale factor in both sectors, though with the 
aforementioned differences in the values. The situation whit x > Xeq is resumed below: 

{const. a < a'^^ , 

a-V2 a'^ <a< a'^^^ , (20) 

a > a'^^, . 

If we recall that the matter-radiation equality for a single sector (ordinary) Universe, 
(fleq)ord) is always bigger than that for a two sectors (mirror) one, (aeq)mir, according to 

(1 + x^) 

(o-eq)!!!!! — ^ ('^eq)ord ^ (Q'eq)ord ; (^1) 

together with our hypothesis x < 1 (from the BBN bound) and (3 > 1 (cosmologically 
interesting situation, i.e., significant mirror baryonic contribution to the dark matter), we 
obtain the useful inequality always verified in a Universe made of O and M sectors 

«b7 < Oeq < ■ (22) 

It's very important to remark (see ref. |2Bf) that at decoupling drops from (p-y/pt) 
to (pb/Pb)- Since p^ oc UyT while pb oc n^T with {n^/rib) ^ 10^ » 1 in the ordinary sector. 
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Figure 3: The trends of the imrror (w^) and ordmary (vg) sound speed compared with the velocity 
dispersion of a typical non baryonic cold dark matter candidate of mass ~ 1 GeV (vcdm)', anr and 
fld indicate the scale factors at which the dark matter particles become non relativistic or decouple. 
The ordinary and mirror models are the same as in figure but the horizontal scale is expanded 
by some decade in order to show the CDM velocity. 

this is a large drop in Vs and consequently in Aj. More precisely, v'^ drops from the value 
{l/3){p^/pb) = {l/3){n^/nb){l + Zdec) to the value (5/3)(rrfecM) = (5/3)(ro/mfe)(l + 
Zdec), with a reduction factor 

Fi^Qbh'^ > 0.026) = = 6.63 • W'^inhh^) , (23) 

(^s)dcc 

where {v^)^^J. and (?^s)dcc indicate the sound speed respectively just before and after 
decoupling. If we consider now the mirror sector and the drop in {v'^)'^ at decoupling, we 
find for the reduction factor 

Fi{x>x,^)=Px-^Fi. (24) 

In the case Ofc/i^ < 0.026, v"^ drops directly from (1/3) to {5/3){Tdec/mb) = (5/3)(ro/mfe)(l+ 
Zrfec) with a suppression 

F2(nbh^ < 0.026) = ^^2^ = 1.9 • 10"^ . (25) 

V^s Jdoc 

It's easy to find that in the mirror sector the reduction factor in the case o^^ > a'^^^ is the 
same as in the ordinary one 

F^{X < Xeq) = F2 . (26) 

Some example: for x = 0.7, F[ 2.9/?Fi; for x = 0.5, F[ = 8/3Fi; for x = 0.3, F{ « 37/3Fi. 
We remark that, if /? > 1, F{ is at least about an order of magnitude larger than Fi. In fact, 
after decoupling {v'J'^ = {5/3){T^^^/mb) = (5/3)(rdec/m6) = (^s)^ (since T^^^ = T^ec), and 
between equipartition and recombination (Ug)^ < (t's)^- The relation above means that 
the drop is smaller in the mirror sector than in the ordinary one. Obviously, before 
equipartition (vi,)'^ = (fg)^ = 1/3, and for this reason the parameter F2 is the same in 
both sectors. 
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In figure|21we plot the trends with scale factor of the mirror sound speed, in comparison 
with the ordinary one. The ordinary model is a typical one with i^^h'^ > 0.026, while the 
mirror model has x = 0.6 and (5 = 2 (this means that M baryonic density is twice the 
ordinary one, chosen in these models about four times its current estimation in order to 
better show the general behaviour). In the same figure we show also the aforementioned 
relative positions of the key epochs (photon-baryon equipartition and decoupling) for both 
sectors, together with the matter-radiation equality. If we reduce the value Q.ifh? , Ob^ 
goes toward higher values, while adoc remains fixed, so that for n^,h? < 0.026 decoupling 
happens before equipartition and the intermediate regime, where Vg oc a~^/^, disappears. 
In an analogous way, if we reduce x, a'^^^ shifts to lower values, until for x < Xeq it occurs 
before the mirror equipartition a'-^^, so that the intermediate regime for vi, disappears. 

In figure 131 the same ordinary and mirror sound speeds are plotted together with the 
velocity dispersion of a typical non baryonic cold dark matter candidate of mass ~ IGeV. 
Note that the horizontal scale is expanded by some decade compared to the figure El 
because the key epochs for the CDM velocity evolution (the epochs when the particles 
become non relativistic, a^r, and when they decouple, Od) occur at much lower scale 
factors. 
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Figure 4: The trends of the mirror Jeans length (Aj) as a function of the scale factor, compared 
with the ordinary Jeans length (Aj). The ordinary model has fl^h^ = 0.08, while the mirror model 
has X = 0.6 and f3 = 2. The horizontal scale is the same as in figure 12 We remark that the same 
behaviours of the ordinary sector are present in the mirror sector for different intervals of scale 
factor. 

3.2 Evolution of the Jeans length and the Jeans mass 

Recalling the definitions Ullf) and (fT^ , and using the results of § 13.11 relative to the sound 
speed, we can compute the evolution of the mirror Jeans length and mass and compare 
them with the analogous quantities for ordinary baryons and CDM. 

We find for the evolution of the adiabatic Jeans length and mass of ordinary baryons 
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in the case > 0.026 



Aj oc 



a 
a 

ai/2 



Mj oc 



oc 



a 

const. 



a < acq 

< a < Odec 

fldcc < « 



(27) 



Otherwise, if Ofe/i^ < 0.026, Ob^ > fldecj there is no intermediate phase < ci < c^dec- 

In the mirror sector it's no more sufficient to interchange Ob-y with a'^^ and Odec with 
a^gj., as made for the sound speed, because from relation H22|) we know that in the mirror 
sector the photon-baryon equipartition happens before the matter-radiation equahty (due 
to the fact that we are considering a mirror sector with more baryons and less photons than 
the ordinary one). It follows that, due to the shifts of the key epochs, the intervals of scale 
factor for the various trends are different. As usual, there are two different possibilities, 
X > Xeq and x < Xeq (which correspond roughly to Qbh'^ > 0.026 and $7^/1^ < 0.026 in 
an ordinary Universe), where, as discussed in S 13.11 for the second one the intermediate 
situationis absent. 

Using the results of § 13.11 for the sound speed, we find the evolution of the adiabatic 
Jeans length and mass in the case x > Xeq 



A't oc 



a 
a 

ai/2 



oc 



a3/2 

const. 

a-3/2 



a < a^^ 

«eq < a < fldec 



^dec 



< a 



(28) 



We plot in figures|l]and[niwith the same horizontal scale the trends of the mirror Jeans 
length and mass compared with those for the ordinary sector; the parameters of both 
mirror and ordinary models are the ones previously used, i.e. Qbh'^ = 0.08, x = 0.6 > x, 
and /? = 2. 



eq 



M, 




Figure 5: The trends of the mirror Jeans mass (Afj) as a function of the scale factor, compared 
with the ordinary Jeans mass (Afj). The models and the horizontal scale are the same as in figure 
il 
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In the ordinary sector the greatest value of the Jeans mass is just before decouphng 
(see ref. ^1), in the interval < a < Cdcc i where 

Mj(a ^ adec) = 1.47 • W'^Mq (1 + (3)'^/^ (n^h^) , (29) 



that for an hypothetical Qh ~ O.l/i ^ and /? = is ~ lO^^M©. Just after decoupling we 
have 

Mj(a > adec) = 2.5 • lO^M© (1 + py^^^ (n^h^) , (30) 

that for Qi, ~ O.lh"'^ and /? = is ~ 10^ Mq. This drop is very sudden and large, changing 
the Jeans mass by F^^^ ~ 1.7 • IQ-^^ {nbh^ f/'^ . 

Otherwise, if Ofe/i^ < 0.026, ab-y > fldeo there is no intermediate phase a^y < a < Odeo 
and Mj(a < Odcc) is larger 

Mj(a < adec) - 3.1 • lO^^^© (1 + /3)-'/' (j^b/^')''^' , (31) 

while after decoupling it takes the value in eq. (|3()|). so that the drop is larger, F2^'^ ~ 8.3 • 
10"^"^. We observe that, with the assumptions = 1 (totally excluded by observations) 
and /? = 0, Mj, 

max (which is the first scale to become gravitationally unstable and collapse 
soon after decoupling) has the size of a supercluster of galaxies. 
If we now consider the expression (|19|) . we have 



<^_/l + /3\/ x4 



/3 y 1 1 + x4 / ' 



(32) 



which can be used to express the value of the mirror Jeans mass in the interval Ceq < 
a < a^gj, (where Mj takes the maximum value) in terms of the ordinary Jeans mass in the 
corresponding ordinary interval ab-y < a < adec- We obtain 



4 \ 3/2 

X 



M'Aa < a'^,,) « (^^—^j . M,{a < adec) , (33) 

which, for P > 1 and x < 1, means that the Jeans mass for the M baryons is lower than 
for the O ones over the entire permitted {P-x) parameter space, with implications for the 
structure formation process. If, e.g., x = 0.6 and /? = 2, then Mj ~ 0.03 Mj. We can also 
express the same quantity in terms oi 0,b, x and /3, in the case that all the dark matter is 
in the form of mirror baryons, as 

/ 4 \ 3/2 

M'jia < a'^J « 3.2 • W'^M^ p-^^l + py^' ( ^ j (n^h^)-^ . (34) 

If we remember eq. (jUJ), we obtain that for the mirror model the drop in the Jeans mass 
at decoupling is (F{)3/2 = p^/^x-^''^{Fif/'^ , which, given our bounds on x and /?, is greater 
than (Fi)3/2. We give here some numerical example: for x = 0.7, (F{)3/2 5/3^/^(Fi)3/2; 
for X = 0.6 and /? = 2 (the case of figures H and , (F{)3/2 « 28(Fi)3/2. 

It's important to stress that these quantities are strongly dependent on the values of 
the free parameters x and /?, which shift the key epochs and change their relative positions. 
We can describe some case useful to understand the general behaviour, but if we want 
an accurate solution of a particular model, we must unambiguously identify the different 
regimes and solve in detail the appropriate equations. 
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Figure 6: The trends of the mirror (Aj) and ordinary (Aj) Jeans length compared with that of a 
typical non baryonic cold dark matter candidate of mass ~ 1 GeV (Appf^^); Onr and ad indicate the 
scale factors at which the dark matter particles become non relativistic or decouple, respectively. 
The models are the same as in figure 01 but the horizontal scale is expanded by some decade to 
show the CDM Jeans length, as in figure 01 



In figures (HI and [7| we plot the trends of the mirror and ordinary Jeans length and 
mass compared with those of a typical non baryonic cold dark matter candidate of mass 
~ IGeV. Apart from the usual expansion of the horizontal scale, due to the much lower 
values of the CDM key epochs as compared to the baryonic ones, a comparison of the 
mirror scenario with the cold dark matter one shows that the maximal value of the CDM 
Jeans mass is several orders of magnitude lower than that for mirror baryons. This implies 
that a very large range of mass scales, which in a mirror baryonic scenario oscillate before 
decoupling, in a cold dark matter scenario would grow unperturbed during all the time 
(for more details see § ]^ . 

For the case x < Xgq, both a'^^ and a'^^^ are smaller than the previous case x > Xeq, 
while the matter-radiation equality remains practically the same; as explained in § (21 the 
mirror decoupling (with the related drop in the associated quantities) happens before the 
matter-radiation equality, and the trends of the mirror Jeans length and mass are the 
following 

( a < a'^^ 

X', OC <! OC OC <^ < « < "dec (35) 

' a a-^ \ const. a^^, < a < 



In this case we obtain for the highest value of the Jeans mass just before decoupling 
the expression 

Mj(a ^a^iJ^ 3.2. 10i^Mer^/2(l + /?)-3/2f^j {Q^h')-^ . (36) 

In the expressions H34|) and (|36() . respectively valid for x > Xeq and x < Xeq, 

are coincident, as we expect. If we consider the differences between the highest mirror 



14 




Figure 7: The trends of the mirror (Afj) and ordmary (Mj) Jeans mass compared with those of a 
typical non baryonic cold dark matter candidate of mass ^ 1 GeV (Mcdm); Onr and Od indicate the 
scale factors at which the dark matter particles become non relativistic or decouple, respectively. 
The model parameters and the horizontal scale are the same as in figure 



Jeans mass for the particular values x = Xeq/2, x = Xeq and x = 2xeq, we obtain the 
following relations: 

^J,max(^cq/2) « (^-j f ^ '^^ 4 J M,max(^cq) « 0.005 Mj^^^JXeq) , (37) 



l + xt 



3/2 



Mj,„,ax(2Xeq) « 2^ ( YT(2^ ) ^J- — ^^^^^ ^ ^J,max(^eq) ■ (38) 

In figure ISl we plot the mirror Jeans mass for the three different possibilities: x < Xeq, 
X > Xeq and x = Xeq (the transition between the two regimes), keeping constant all other 
parameters. In these three models the matter-radiation equality is the only key epoch 
which remains almost constant. The change in the trends when x becomes lower than 
Xeq, due to the fact that a'^^^ becomes lower than Oeq, generates an evident decrease of the 
Jeans mass. 

3.3 Evolution of the Hubble mass 

The trends of the Hubble length and mass are expressed, as usually, by 



a 



2 ( „3 



(Ah) \ a < aeq 



Ah oc < 0/2 Mh oc — ^ Of ^ s/2 (39) 



"eq 



It should be emphasized that, as for the ordinary baryons, during the period of dom- 
ination of photons (a < a{,^) the mirror baryonic Jeans mass is of the same order of the 
Hubble mass. In fact, following definitions we find 

M=W)!.26. (40) 

Mh A^ 

We plot the trends of the Hubble mass in figures and together with other fun- 
damental mass scales. 
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Figure 8: The trends of the mirror Jeans mass for the cases x < x^q (solid line), x — Xcq (dotted) 
and X > Xcq (dashed). The model with x > Xcq is the same as in figure El the others are obtained 
changing only the value of x and keeping constant all the other parameters. As clearly shown in 
the figure, the only key epoch which remains almost constant in the three models is the matter- 
radiation equality; the mirror baryon-photon equipartition and decoupling indicated are relative 
to the model with x < Xcq. It's also evident the change in the trends when x becomes lower than 
Xcq, due to the fact that a^^^ becomes lower than Ooq. 

3.4 Dissipative effects: collisional damping 

A peculiar feature of the mirror baryonic scenario is that mirror baryons undergo the colli- 
sional damping as ordinary ones. This dissipative process modify the purely gravitational 
evolution of perturbations. The physical phenomenon is the interaction between baryons 
and photons before the recombination, and the consequent dissipation due to viscosity 
and heat conduction. Around the time of recombination the perfect fluid approximation 
breaks down, and the perturbations in the photon-baryon plasma suffer from collisional 
damping. As decoupling is approached, the photon mean free path increases and photons 
can diffuse from the overdense into the underdense regions, thereby smoothing out any 
inhomogeneities in the photon-baryon plasma. This effect is known as Silk damping |3nj . 

In order to obtain an estimate of the effect, we follow ref. jSIj for ordinary baryons, 
and then we extend to mirror baryons. We consider the photon mean free path 



where is the electron ionization factor, oc a is the number density of the free 
electrons and ctt is the cross section for Thomson scattering. Clearly, photon free stream- 
ing should completely damp all baryonic perturbations with wavelengths smaller than A-y. 
Damping, however, occurs on scales much larger than since the photons slowly diffuse 
from the overdense into the underdense regions, dragging along the still tightly coupled 
baryons. Integrating up to decoupling time we obtain the total distance traveled by a 
typical photon 




(41) 



As 




) 



3/4 



Mpc 



(42) 
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and the associated mass scale, the Silk mass, given by 



4 /As 
Ms = gTTp, ( - 



(43) 



which, assuming O^/i^ ~ 0.02, gives Ms ~ 8 x 10-^^ M0. This dissipative process causes 
that fluctuations on scales below the Silk mass are completely washed out at the time of 
recombination and no structure can form on these scales. This has consequences on large 
scale structure power spectrum, where small scales have very little power. 
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Figure 9: The trend of the mirror Silk mass (Mg) over a cosmologically interesting range of x, 
which contains Xeq (wc considered Vt^h^ ~ 0.15, so Xeq — 0.3). The axis are both logarithmic. We 
show for comparison also the values of the ordinary Silk mass (Ms) and of the free streaming mass 
(Mfs) for typical HDM and WDM candidates. 



In the mirror sector too, obviously, the photon diffusion from the overdcnsc to under- 
dense regions induces a dragging of charged particles and washes out the perturbations at 
scales smaller than the mirror Silk scale 



A^~3/(a;)(/3J^b/i2)-3/4 Mpc , 



(44) 



where f{x) = x^/^ for x > Xeq and f{x) = {x/xec^^^'^x^ for x < Xeq. Thus, the density 
perturbation scales running the linear growth after the matter-radiation equality epoch 
are limited by the length Ag. The smallest perturbations that survive the Silk damping 
will have the mass 

Mk ~ [/(x)/2]3(/3J^,/,2)-5/4ioi2 ^ (45) 

which should be less than 2 x 10^^ Mq in view of the BBN bound x < 0.64. Interestingly, 
for X ^ a:;eq we obtain, for the current estimate of i^mh^ and if all the dark matter is made 
of mirror baryons, M'g ~ 10"'^'' Mq, a typical galaxy mass. 

At this point it is very interesting a comparison between different damping scales, 
collisional (ordinary and mirror baryons) and collisionless (non-baryonic dark matter) . We 
know that for hot dark matter (as a neutrino with mass ~10 eV) Mpg 10^^ Mq, while 
for a typical warm dark matter candidate with mass ~1 keV, M^^^ ~ 10^ — lO^'^ Mq. 
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From eq. H45() it is evident that the dissipative scale for mirror Silk damping is analogous 
to that for WDM free streaming. Consequently, the cutoff effects on the corresponding 
large scale structure power spectra are similar, though with important differences due to 
the presence of oscillatory features, which makes them distinguishable one from the other 
(for a detailed presentation of these power spectra see the Paper 2 In figure |H1 we 

show this comparison together with the trend of the mirror Silk mass over a cosmologically 
interesting range of x. 

4 Scenarios 

After the description of the fundamental scales for structure formation, let us now collect 
all the informations and discuss the mirror scenarios. They are essentially two, according 
to the value of x, which can be higher or lower than Xeq, and are shown respectively in 
figures Uni and ITTl which will be our references during the present section. 

Typically, adiabatic perturbations for mirror baryons with sizes larger than the max- 
imum value of the Jeans mass, which is Mj(aeq) for x > Xeq and Mj(a^g^) for x < Xeq, 
experience uninterrupted growth. In particular, they grow as oc before matter- 
radiation equality and as (5f, oc o after equality. Fluctuations on scales in the mass interval 
,max grow as 5ft OC (2 while they are still outside the Hubble radius. After 
entering the horizon and until recombination these modes oscillate like acoustic waves. 
The amplitude of the oscillation is constant before equilibrium but decreases as a~^^^ 
between equipartition and recombination. After decoupling the modes become unstable 
again and grow as Sb oc a. Finally all perturbations on scales smaller than the value of the 
Silk mass are dissipated by photon diffusion. 




Figure 10: Typical evolution of a perturbed scale Mport in adiabatic mirror baryonic dark matter 
scenario with x > Xcq- The figure shows the Jeans mass Mj, the Silk mass Mg and the Hubble mass 
Mh- The time of horizon crossing of the perturbation is indicated by flent- Are also indicated the 
three evolutionary stages: during stage A (a < Ocnt < Ocq) the mode grows as Sb o: a^; throughout 
stage B (oent < a < a'^^^) the perturbation oscillates; finally, in stage C (a > a'^^^) mode 
becomes unstable again and grows as oc a. Fluctuations with size smaller than Mg are wiped 
out by photon diffusion. 



18 



Given this general behaviour, the schematic evolution of an adiabatic mode with a 
reference mass scale Mpert, with Mg < Mpert < Afj(aeq), is depicted in figure ITUl for 
X > Xgq. We distinguish between three evolutionary stages, called A, B and C, depending 
on the size of the perturbation and on the cosmological parameters il^h'^, x and /3, which 
determine the behaviour of the mass scales, and in particular the key moments (time of 
horizon crossing and decoupling) and the dissipative Silk scale. During stage A, i.e. before 
the horizon crossing (a < acnt < «eq)) the mode grows as 6b oc a^; throughout stage B 
(«ent < a < «dec) perturbation enters the horizon, baryons and photons feel each 
other, and it oscillates; finally, in stage C (a > a'^cc)^ the photons and baryons decouple, 
and the mode becomes unstable again growing as 6b oc a. We remark that fluctuations 
with sizes greater than Mj(aeq) grow uninterruptedly (because after horizon crossing the 
photon pressure cannot balance the gravity), changing the trend from before MRE to 
a after it, while those with sizes smaller than Mg are completely washed out by photon 
diffusion. 

After decoupling, all surviving perturbations (those with Mpert > -^s) grow steadily 
until their amplitude becomes so large that the linear theory breaks down and one needs 
to employ a different type of analysis. 




Figure 11: Typical evolution of a perturbed scale Mpert in adiabatic mirror baryonic dark matter 
scenario with x < Xoq- The value of Mpert is the same as in figure 11171 The time of horizon crossing 
of the perturbation is indicated by Oent- The figure shows the Jeans mass Mj, the Silk mass Mg 
and the Hubble mass Mh. Unlike the case x > Xcq (shown in the previous figure), now there are 
only the two evolutionary stages A (a < aent) and C (a > flcnt). Fluctuations with size smaller 
than Mg are wiped out by photon diffusion, but in this case the Silk mass is near to the maximum 
Jeans mass. 

If we look, instead, at the schematic evolution of an adiabatic mode with the same 
reference mass scale Mpert but for x < x^q, as reported in figure 1111 we immediately 
notice the lower values of the maximum Jeans mass and the Silk mass, which are similar. 
Therefore, for the plotted perturbative scale there are now only the two stages A and 
C. In general, depending on its size, the perturbation mass can be higher or lower than 
the Silk mass (and approximatively also than the maximum Jeans mass), so modes with 
-^pert > -^s grow continuously before and after their horizon entry, while modes with 
-^pert < -^s completely washed out. 
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We find that Mj becomes smaller than the Hubble horizon mass Mh starting from a 
redshift 

Zc = 3750 x-^ {^mh'^) , (46) 

which is about Zeq for x = 0.64, but it sharply increases for smaller values of as shown in 
figure n We can recognize this behaviour also watching at the intersections of the lines for 
Mj and Mh in figures ITUl and ITTl Thus, density perturbation scales which enter horizon 
at 2: ~ Zeq have masses larger than M'j and thus undergo uninterrupted linear growth 
immediately after teq- Smaller scales for which Mj > Mjj would instead first oscillate. 
Therefore, the large scale structure formation is not delayed even if the mirror decoupling 
did not occur yet, i.e. even if x > Xcq. 

When compared with non baryonic dark matter scenarios, the main feature of the mir- 
ror baryonic scenario is that the M baryon density fluctuations should undergo the strong 
coUisional damping around the time of M recombination, which washes out the pertur- 
bations at scales smaller than the mirror Silk scale. It follows that density perturbation 
scales which undergo the linear growth after the MRE epoch are limited by the length 
Ag. This could help in avoiding the excess of small scales (of few Mpc) in the CDM power 
spectrum without tilting the spectral index. To some extent, the cutoff effect is analogous 
to the free streaming damping in the case of warm dark matter (WDM), but there are 
important differences. The point is that, alike usual baryons, the MBDM shows acoustic 
oscillations with an impact on the large scale structure (LSS) power spectrum [Uljnj. In 
particular, it is tempting to imagine that the M baryon oscillation effects are related to 
the anomalous features observed in LSS power spectra. 

In addition, the MBDM oscillations transmitted via gravity to the ordinary baryons, 
could cause observable anomalies in the CMB angular power spectrum for Ts larger than 
200. This effect can be observed only if the M baryon Jeans scale Aj is larger than the 
Silk scale of ordinary baryons, which sets a principal cutoff for CMB oscillations around 
I ~ 1200. This would require enough large values of x, near the upper bound fixed by the 
BBN constraints, and, together with the possible effects on the large scale power spectrum, 
it can provide a direct test for the MBDM (verifiable by the higher sensitivity of next CMB 
and LSS observations). For a complete discussion of the CMB and LSS power spectra for 
a Mirror Universe see the Paper 2 22> . 

Clearly, for small x the M matter recombines before the MRE moment, and thus it 
behavesas the CDM as far as the large scale structure is concerned. However, there still 
can be crucial differences at smaller scales which already went non-linear, like galaxies. 
In our scenario, dark matter in galaxies and clusters can contain both CDM and MBDM 
components, or can be even constituted entirely by the mirror baryons. 

One can question whether the MBDM distribution in halos can be different from that 
of the CDM. Simulations show that the CDM forms triaxial halos with a density profile 
too clumped toward the center, and overproduces the small substructures within the halo. 
Since MBDM constitutes a kind of collisional dark matter, it may potentially avoid these 
problems, at least the one related with the excess of small substructures. 

It's also worth noting that, throughout the above discussion, we have assumed that the 
matter density of the Universe is close to unity. If, instead, the matter density is small and 
a vacuum density contribution is present, there is an additional complication due to the 
fact that the Universe may become curvature dominated starting from some redshift -^curv • 
Given the current estimate — 0.7, this transition has yet occurred and the growth of 
perturbations has stopped around Zcurv, when the expansion became too rapid for it. 
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5 Evolution of perturbations 



As a result of the studies done in previous sections, in this section we finally consider the 
temporal evolution of perturbations, as function of the scale factor a. All the plots are 
the results of numerical computations obtained using a Fortran code originally written for 
the ordinary Universe and modified to account for the mirror sector. 

We used the synchronous gauge and the evolutionary equations presented in ref. P^ . 
The difference in the use of other gauges is limited to the gauge-dependent behaviour of 
the density fluctuations on scales larger than the horizon. The fluctuations can appear as 
growing modes in one coordinate system and as constant mode in another, that is exactly 
what occurs in the synchronous and the conformal Newtonian gauges. 

In the figures we plot the evolution of the components of a Mirror Universe, namely 
the cold dark matter^, the ordinary baryons and photons, and the mirror baryons and 
photons, changing some parameter to evaluate their influence. 

First of all, we comment figure IT^ . which is the most useful to recognize the general 
features of the evolution of perturbations. Starting from the smallest scale factor, we 
see that all three matter components and the two radiative components grow with the 
same trend (as a^), but the radiative ones have a slightly higher density contrast (with a 
constant rate until they are tightly coupled); this is simply the consequence of considering 
adiabatic perturbations, which are linked in their matter and radiation components by the 
adiabatic condition (|lfl|) . This is the situation when the perturbation is out of horizon, 
but, when it crosses the horizon, around a ~ 10~^, things drastically change. Baryons 
and photons, in each sector separately, become causally connected, feel each other, and 
begin to oscillate for the competitive effects of gravity and pressure. Meanwhile, the CDM 
density perturbation continues to grow uninterruptedly, at first reducing his rate from 
to In a (due to the rapid expansion during the radiation era), and later, as soon as MRE 
occurs (at a ~ 3 X 10^^ for the considered model), increasing proportionally to a. The 
oscillations of baryons and photons continue until their decoupling, which in the mirror 
sector occurs before than in the ordinary one (scaled by the factor x). This moment is 
marked in the plot as the point where the lines for the two components move away one from 
the other. From this point, the photons in both sectors continue the oscillations until they 
are completely damped, while the M and O baryons rapidly fall into the potential wells 
created by the cold dark matter and start growing as a. We remark that it's important the 
way in which the oscillation reaches the decoupling; if it is compressing, first it slows down 
(as if it continues to oscillate, but disconnected from the photons), and then it expands 
driven by the gravity; if, otherwise, it is expanding, it directly continues the expansion 
and we see in the plot that it immediately stops to oscillate. In this figure we have the 
first behaviour in the mirror sector, the second one in the ordinary sector. 

In figures I12| 1131 and we compare the behaviours of different scales for the same 
model. The scales are given by log A;(Mpc~^) = —0.5, —1.0, —1.5, —2.0, —2.5, —3.0, where 
= 27r/A is the wave number. The effect of early entering the horizon of small scales (those 
with higher k) is evident. Going toward bigger scales the super horizon growth continue 
for a longer time, delaying more and more the beginning of acoustic oscillations, until it 
occurs out of the coordinate box for the bigger plotted scale (log A; = —3.0). Starting 
from the scale log/c = —1.5, the mirror decoupling occurs before the horizon entry of 
the perturbations, and the evolution of the mirror baryons density is similar to that of 
the CDM. The same happens to the ordinary baryons too, but for logfc ^ —2.0 (since 
they decouple later), while the evolution of mirror baryons is yet indistinguishable from 

^As non baryonic dark matter we consider only the cold dark matter, which is at present the standard 
choice in cosmology. 
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that of the CDM. For the bigger scales (log A; < —2.5) the evolution of all three matter 
components is identical. 

As previously seen, the decoupling is a crucial point for structure formation, and it 
assumes a fundamental role specially in the mirror sector, where it occurs before than in 
the ordinary one: mirror baryons can start before growing perturbations in their density 
distribution. For this reason it's important to analyze the effect of changing the mirror 
decoupling time, obtained changing the value of x and leaving unchanged all other param- 
eters, as it is possible to do using figures IT^. [T31 and [TBI for x = 0.6,0.5,0.4,0.3,0.2 and 
the same scale log k = —1.0. It is evident the shift of the mirror decoupling toward lower 
values of a when reducing x, according to the law Q, which states a direct proportionality 
between the two. In particular, for x < x^q ~ 0.3 mirror decoupling occurs before the 
horizon crossing of the perturbation, and mirror baryons mimic more and more the CDM, 
so that for x ~ 0.2 the perturbations in the two components are indistinguishable. For 
the ordinary sector apparently there are no changes, but at a more careful inspection we 
observe some difference due to the different amount of relativistic mirror species (propor- 
tional to x^), which slightly shifts the matter-radiation equality. This effect is more clear 
in figure ITTl where we plot only the CDM and the ordinary baryons for the cases x = 0.2 
and 0.6: for the lower value of x there are less mirror photons, the MRE occurs before 
and the perturbation in the collisionless component starts growing before proportionally 
to the scale factor; thus, when the baryons decouple, their perturbation rapidly grows to 
equalize that in the CDM, which meanwhile has raised more for the lower x. 

Obviously, these are cases where the CDM continues to be the dominant form of dark 
matter, and drives the growth of perturbations, given its continuous increase. In any case, 
if the dominant form of dark matter is made of mirror baryons the situation is practically 
the same, as visible comparing figures H^b and II 8b . (where we see only slight differences on 
the CDM and mirror baryons behaviours in the central region of the plots), since mirror 
baryons decouple before than ordinary ones and fall into the potential wells of the CDM, 
reinforcing them. 

Finally, in the interesting case where mirror baryons constitute all the dark matter, 
they drive the evolution of perturbations. In fact, in figure ITHb we clearly see that the 
density fiuctuations start growing in the mirror matter and the visible baryons are involved 
later, after being recombined, when they rewrite the spectrum of already developed mirror 
structures. This is another effect of a mirror decoupling occurring earlier than the ordinary 
one: the mirror matter can drive the growth of perturbations in ordinary matter and 
provide the rapid growth soon after recombination necessary to take into account of the 
evolved structures that we see today. 

After all the considerations made in this paper, it is evident that the case of mirror 
baryons is very interesting for structure formation, because they are collisional between 
themselves but collisionless for the ordinary sector, or, in other words, they are self- 
collisional. In this situation baryons and photons in the mirror sector are tightly coupled 
until decoupling, and structures cannot grow before this time, but the mirror decoupling 
happens before the ordinary one, thus structures have enough time to grow according 
to the limits imposed by CMB and LSS (something not possible in a purely ordinary 
baryonic scenario). Another important feature of the mirror dark matter scenario is that, 
if we consider small values of x, the evolution of primordial perturbations is very similar 
to the CDM case, but with a fundamental difference: there exist a cutoff scale due to the 
mirror Silk damping, which kills the small scales, overcoming the problems of the CDM 
scenario with the excessive number of small satellites. These are important motivations 
to go further in this work and investigate the effects of the mirror sector on the CMB and 
LSS power spectra, as we will do in Paper 2 [23] • 
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Figure 12: Evolution of perturbations for the components of a Mirror Universe: cold dark matter 
(solid line) , ordinary baryons and photons (dot-dashed and dotted) and mirror baryons and photons 
(long dashed and dashed). The model is a flat Universe with = 0.3, fib/i^ = 0.02, = 0.04 
(/3 = 2), /i = 0.7, X = 0.6, and plotted scales are logfc(Mpc~^) = -0.5 (a) and -1.0 (6). 
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^ure 13: The same as in figure IT^ but for scales logfc(Mpc^"^) = —1.5 (a) and —2.0 
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Figure 14: The same as in figure IT^ but for scales logfc(Mpc ^) — —2.5 (a) and —3.0 
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Figure 15: Evolution of perturbations for the components of a Mirror Universe: cold dark matter 
(solid line), ordinary baryons and photons (dot-dashed and dotted) and mirror baryons and photons 
(long dashed and dashed). The model is a flat Universe with n,„ = 0.3. ftbh'^ = 0.02, = 0.04 
[P = 2),h = 0.7, X = 0.5 (a) or 0.4 (6), and plotted scale is logA;(Mpc"^) = -1.0. 
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Figure 16: The same as in figure El but for x ~ 0.3 (a) and 0.2 (6). 
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Figure 17: Evohition of perturbations in a Mirror Universe for cold dark matter (solid and long 
dashed lines) and ordinary baryons (dashed and dotted). The models are flat, with = 0.3, 
fife/i^ = 0.02, = 0.04 (/? = 2), /i = 0.7, x = 0.2 (solid and dashed) and 0.6 (long dashed and 
dotted), and logfc(Mpc~^) = —1.0. The highest scale factor plotted here is 0.05, while in all other 
flgures it is 0.02. 
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Figure 18: Evolution of perturbations for the components of a Mirror Universe: cold dark matter 
(solid line), ordinary baryons and photons (dot-dashed and dotted) and mirror baryons and photons 
(long dashed and dashed). The models arc flat with = 0.3, flbh^ = 0.02, = 40;, (a) or 
(n„ - (no CDM) (6), h = 0.7, x = 0.4, and logA;(Mpc-^) = -1.0. 
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6 Conclusions 



We have studied the structure formation in Hnear regime for a Universe with mirror dark 
matter. This topic was introduced in previous works I', and here we analyzed the details 
of the evolution of all the relevant scales (sound speed, Jeans length and mass, dissipative 
scales) comparing with the ordinary and the CDM cases and for different mirror scenarios, 
together with the evolution of perturbations varying a lot of parameters. 

We extended the Jeans theory to the more general case where the Universe is made of 
two sectors, ordinary and mirror, using a pre-existent bound (coming from the BBN limits) 
on the mirror temperature: in the relativistic expansion epoch the cosmological energy 
density is dominated by the ordinary component and the mirror one gives a negligible 
contribution, while for the non-relativistic epoch the complementary situation can occur 
when the mirror baryon matter density is bigger than the ordinary one, fi^ > r^f,, and 
hence the mirror baryonic dark matter (MBDM) can contribute the dark matter of the 
Universe along with the CDM or even constitute a dominant dark matter component. 

Since the existence of a mirror hidden sector changes the time of key epochs, there 
are important consequences on the structure formation scenario for a Mirror Universe. 
We studied this scenario in presence of adiabatic scalar density perturbations, which are 
at present the most probable kind of primordial fluctuations, in the context of the Jeans 
gravitational instability theory. 

Given that the physics is the same in both sectors, key differences are the shifts of 
fundamental epochs, namely the baryon-photon equipartition and the matter-radiation 
decoupling occur in the mirror sector before than in the ordianary one: a'-^^ < ab-^; a'^^^ < 
Odec- The first step was the study of the mirror sound speed and its comparison with 
the ordinary one and with the velocity dispersion of a typical cold dark matter candidate. 
From this study we obtained the mirror Jeans length and mass, again to be compared with 
the same quantities obtained for the ordinary sector and for the cold dark matter. There 
are two different possibilities, according to the value of x, which can be higher or lower than 
Xeq ~ 0.046(ilm/i^)~^, the value for which mirror decoupling occurs at matter-radiation 
equality time. 

The values of the length and mass scales clearly depend on the mirror sector tem- 
perature and baryonic density, but we found that M'j is always smaller than Mj, with a 
typical ratio ~ 50 for x > Xeq, while for cold dark matter it is several orders of magnitude 
lower. For x < x^.^, M'j is few orders smaller than the same quantity obtained for the case 

X > ^eq- 

Another important quantity to describe the structure evolution is the dissipative scale, 
represented in the mirror sector by the mirror Silk scale. We found that it is much lower 
than the ordinary one, obtaining Mg ~ 10^*^ Mq (around a typical galactic mass) for 
X ~ Xeq, a value similar to the free streaming scale for a typical warm dark matter 
candidate, and much higher than the one for cold dark matter. 

We collected all these informations in two different mirror scenarios, for x > Xeq and 
X < Xeq. In the latter case we obtained a maximum mirror Jeans scale similar to the 
mirror Silk mass, so that practically all perturbations with masses greater than this Silk 
mass grow uninterruptedly, just as in a cold dark matter scenario. 

After this, we modified a numerical code existing for the standard Universe in order to 
take into account a hidden mirror sector, and computed the evolution of perturbations in 
the linear regime for the components present in a Mirror Universe, namely the ordinary 
and mirror baryons and photons, and the cold dark matter. We did this for various mirror 
temperatures and baryon densities, and for different perturbative scales, finding all the 
features predicted by our structure formation study (as for example the mirror decoupling 
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and the CDM-like behaviour for low values). 

This was the first step in order to obtain the photon and matter power spectra for a 
Universe with mirror dark matter and compare them with observations, and this will be 
exactly the object of the next paper of this series. 
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